hep-ph 9403361 

SLAC-PUB-6464 
Revised Version 
May 1994 



ANGULAR DISTRIBUTIONS IN THE DRELL-YAN PROCESS: 
A CLOSER LOOK AT HIGHER TWIST EFFECTS 



A. Brandenburg 1,2 , S. J. Brodsky 2 , V. V. Khoze 2 and D. Miiller 2 ' 3 

Stanford Linear Accelerator Center 
Stanford University, Stanford, California 94309 



Abstract 

We calculate the angular distribution of the lepton produced in the Drell-Yan reaction 
taking into account pion bound state effects. We work in the kinematic region where one 
of the pion constituents goes far off-shell, which allows us to treat the bound state problem 
perturbatively. We show that the angular distribution is very sensitive to the shape of the 
pion distribution amplitude. The model we discuss fits the data if we choose a two-humped 
pion distribution amplitude suggested by QCD sum rules. 



Submitted to Phys. Rev. Let. 



1 Max Kade fellow 

2 Work supported by the Department of Energy, contract DE-AC03-76SF00515 

3 Supported by Deutscher Akademischer Austauschdienst 



Lepton pair production in hadron-hadron collisions provides a basic testing ground 
for our understanding of strong interactions. Extensive experimental and theoretical work 
has been done in the past two decades (for reviews see Ref. [1]). In particular, the angular 
distribution of the lepton pair has been studied in detail, revealing a fatal disagreement 
of the QCD improved parton model prediction [2,3] with the data [4-6]. Recently it has 
been proposed that this problem may be resolved if the nontrivial structure of the QCD 
vacuum induces spin correlations between the initial state partons [7]. In this letter we 
pursue another way to go beyond the standard parton model picture, namely we consider 
contributions to the angular distribution induced by hadron bound state effects. Our 
approach is close in spirit to the higher twist model of [8,9]. 

The angular distribution of the fx + in 

7T- +N +X -> fj+ + fT +X (1) 

may be parameterized in general as follows: 

-^r ~ 1 + A cos 2 9 + \i sin 29 cos </> + - sin 2 9 cos 20. (2) 
<j ail 2 

Here 9 and are angles defined in the muon pair rest frame and A, fx and v are angle- 
independent coefficients. The naive parton model (Drell-Yan picture [10]) views the pro- 
duction of the virtual photon 7* in (1) as originating from the annihilation of two un- 
corrected constituent quarks, resulting in an angular distribution of the form 1 + cos 2 9. 
This result follows simply from the fact that the virtual photon is produced transversely 
polarized in the annihilation of two on-shell fermions. 

In order to describe the boson transverse momentum distribution d 2 o jdQ\ one has 
to take into account radiative corrections to the Drell-Yan model. The QT-distribution 
has been calculated in the QCD-improved parton model to the order of 0(a s ) with re- 
summation of the soft gluons at the leading double logarithmic accuracy (see Ref. [11] and 
references therein). This approach was used in [3] to compute the angular distribution at 
fixed transverse momentum. The deviations from the 1 + cos 2 9 behavior were found to 
be less than 5% in the range < Qt < 3 GeV [3]. However, the NA-10 measurements 
from CERN [4] and the Chicago-Iowa-Princeton collaboration [5,6] show a quite different 
behavior. In the limit where the momentum fraction x of one of the pion constituents is 
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very close to 1 and for moderate transverse momenta of the muon pair, the value of A 
turns strongly negative [6] , consistent with a sin 2 9 distribution. This implies that in this 
kinematic limit the virtual photon is produced with longitudinal polarization, rather than 
transverse. Furthermore, the data [4-6] is observed to have a strong azimuthal modulation 
(nonzero \i and v in (2)), an effect which is missing in standard QCD. The Lam-Tung sum 
rule [2], 1 — A — 2u = 0, which follows from the approach used in [3] is also badly violated 
by the experimental data. 

One way to go beyond the standard treatment is to take into account the pion bound 
state effects [8,9]. We want to treat the bound state problem perturbatively; thus we will 
restrict ourselves to a specific kinematic region in which the momentum fraction x of one 
of the pion constituents is large, x > 0.5. In fact, in the large x region the off-shell nature 
of the annihilating quark from the projectile is crucial, and thus the operative subprocess 
must involve the correlated multi-parton structure of the projectile. In effect the dominant 
subprocess in the off-shell domain is n~q — > fi + (i~q. We resolve the pion by a single hard 
gluon exchange [12]. The main contribution to reaction (1) then comes from the diagrams 
of Fig. la,b [8,9]. We see from diagram la that the u quark propagator is far off-shell, 
p? = — Qt/(1 — x u)- The second diagram is required by gauge invariance^. The leading 
contribution to the amplitude M for the reaction 

u + tv- ^7*+X^ / u+ + /U -+X (3) 

is obtained [12] by convoluting the partonic amplitude T(u + ud — > 7* + d — > fx + + yT + d) 
with the pion distribution amplitude 4>(z,Q 2 ), 

M = [ dz (p(z,Q 2 ) T, (4) 
Jo 

where Q 2 ~ <5t/(1 — x ) is the cutoff for the integration over soft momenta in the definition 
of 4>. In the regime where Q 2 and Q 2 are compatible, one cannot use the usual probabilistic 
factorization of the structure functions and the hard annihilation subprocess [13]. 
For the hadronic differential cross section we have 

t In a physical gauge the contribution of the second diagram is purely higher twist, that 
is it contains extra powers of \/Q\jQ 2 - 
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+x u -(5tS _1 (1 -Xu)' 1 ) S(Q 2 - sx u Xu + Qt{1 - Xu) -1 ) + {u — > d,u^ d}. (5) 

Here Q M is the four-momentum of 7* in the hadronic center of mass system, x u ^) is 
the light-cone momentum fraction of the u(u) quark and G u /n is the parton distribution 
function of the nucleon. The longitudinal momentum fraction of the photon is defined as 
xl = 2Ql I y/s and it should be noted that its maximum value, x™ ax = 1 — s~ 1 (Q 2 + 2 Q T ) 
is slightly less than 1. The second term on the right hand side of (5) is the same as the first 
one with quark flavors interchanged. This term gives the contribution from the nucleon 
sea. In Fig. lc we show a typical contribution to the hadronic cross section. 

We note that no primordial or intrinsic transverse momenta have been introduced. 
The single gluon exchange is the only source of Qt in the model discussed. We also 
neglected the quark masses and the mass of the projectile which are small compared to Q. 

In analogy to eq. (2) we parameterize the angular distribution as follows, 

Q 2 da ( Q 2 da 



dQ 2 dQ\dxLdVt \dQ 2 dQ^dxL 

3 1 is 

■(1 + A cos 2 6 + fj, sin26>cos</> + — sin 2 cos 20), (6) 



4ttA + 3 v ^ r 2 

where the angular distribution coefficients A, \i and v are now functions of the kinematic 
variables xl, Qt/Q 2 and Q 2 /s. 

We work in the Gottfried- Jackson frame where the z axis is taken to be the pion 
direction in the muon pair rest frame and the y axis is orthogonal to the ir~N plane. With 
some algebra, using eqs. (4)-(6), we arrive at the following expression for A, ji and v, 

X(x, Qt/Q 2 ) = 2N- 1 j (1 - x) 2 [(1m I(x)) 2 + (F + Re I(x)) 2 ] 

-4Q 2 T /Q 2 x 2 F 2 + Q 4 T /Q 4 x 2 F 2 j, (7) 
fi(x, Q 2 T /Q 2 ) = -4N- 1 ^Q 2 ,/Q 2 Fx!^ (1 — x)[F + Re I(x)] + Q 2 T /Q 2 xF j, (8) 



v(x, Qt/Q 2 ) = SN^Ql/Q 2 ^! - x)F [F + Re I(x)} , (9) 



where 



F= f 1 dz t^l, (10) 

Jo z 

m-f 1 ** , .. . (11) 
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and 



z (z + x — 1 + ie) ' 
N(x, Q 2 T /Q 2 ) = 2<j (1 - £) 2 [(Im 7(£)) 2 + (F + Re /(x)) 2 ] 



+4Q 2 7Q 2 x 2 F 2 + /Q 4~2 F 2 | (12) 

The variable x, acts to resolve the distribution amplitude much like the Bjorken variable 
resolves the structure functions, 

xu _l x L + ^x 2 L +4s- 1 (Q 2 + Q' 2 T ) 
X ~1 + Q 2 T /Q 2 2 l + Q 2 T /Q 2 ' [ ' 

The factors 1/z in eq. (10), (11) come from the gluon propagators and the factors l/(z + 
x — 1 ± ie) arise from the quark propagator of Fig. lb. 

In contrast to Refs. [9] and [14] we did not omit terms 0((1 — Xy)Q\jQ 2 ) and 
0((1 — Xu)' 1 Qt/Q 4 ) an d of higher orders. The nucleon distribution function G q / N does 
not appear in (7)- (9); thus only the pion distribution amplitude <f>(z, Q 2 ) has to be specified. 
Before doing so we give some technical comments on our calculation. 

We note that the internal quark line of Fig. lb can go on-shell. The amplitude M of 
equation (4) , however, is always regular due to the ^-integration [14] for realistic choices of 
(f)(z, Q 2 ). This also can be read off from (11). The fact that the internal line goes on-shell 
does not cause a Sudakov suppression since our diagrams are the lowest order contribution 
of an inclusive process. In other words gluon emission to the final state will occur in 
the higher order corrections. Only when x u approaches unity, where gluon emission is 
prohibited by kinematics, the Sudakov suppression will arise. 

Our model and the parton model are not complementary, but rather different approx- 
imations to the Drell-Yan process. The diagrams of Fig. la,b give the whole leading order 
contribution in the specific kinematic region of large enough x u , x u > 0.5 [9]. This is so 
because the gluon exchange is the resolution of the pion bound state and not a radiative 
correction. 
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Now we can present our final results for A, n and v for different choices of the pion 
distribution amplitude <f>(z, Q 2 ). We find in general that the values of \i and v are very- 
sensitive to the choice of <f>(z, Q 2 ) which we always take to be positive, symmetric, i.e. 
(p(z,Q 2 ) = <p{l — z,Q 2 ), and normalized, f Q dz <p(z,Q 2 ) = 1. Thus we will not restrict 
ourselves to the simplest case of <j)(z) ~ S(z — 1/2) considered in [9]. 

This sensitivity can be illustrated for the special case of x = 0.5 for which Re I = 
-2F. From eqs. (7)-(9) we get, 

\{x - 0.5, Q T /Q + 4yr2fl2 + AQl/Q2 + Q , /Q4 , (14) 



»{x - 0.5, Q T /Q + 4n2a2 + 4Q ^ /Q2 + Q , /Q4 , (15) 

^ - °- 5 ' ^ - i + w + 4 S% + ^ ' (l6) 

where a = <p(z = 0.5) /F. 

From these formulas we see that \x and v in this case are not suppressed only if a is a 
sufficiently small number. This is true for the so-called two-humped distribution amplitude 
[15] which has a dip around z = 0.5. On the other hand, the choice of a convex distribution 
amplitude, e.g. the asymptotic one, <J)(z) = 6z(l — z), will always produce suppressed \x 
and v at x = 0.5. 

We return now to our general results, eqs. (7)-(9). In Fig. 2 we plot A, fx, v and 
2v — (1 — A) versus Xu for a/ Q\IQ 2 = 0.25 for different choices of 4>(z, Q 2 ) together with 
the data of Ref. [5]. For the two- humped distribution amplitude we have chosen the 
evolution parameter Q 2 to be effectively ~ 4 GeV 2 . The solid line is the result for the 
two-humped <f>(z) where powers of {Q\IQ 2 Y i l 2 were dropped for n > 3 in eqs. (7)-(9). We 
note that corrections to our model may induce such terms, thus the difference between the 
dashed and the solid lines should be viewed as the uncertainty of our predictions. We also 
show the data points of Ref. [5] averaged in the intervals 4.05 < ^/Q 2 < 8.55 GeV and 
< y/Q% < 5 GeV. 

In Fig. 3 the same quantities are shown versus \JQ\ for x^ = 0.6 and ^/Q 2 = 6 GeV. 
The data points in this case are averaged over intervals 4.05 < \[Q 2 < 8.55 GeV and 
0.2 < Xu < 1 and taken from Ref. [5]. All the data points were averaged over the intervals 



defined above in Ref. [5]. We would rather prefer to use the unaveraged data which are 
not available. The use of the averaged over x^ data in Fig. 3 forced us to fix the value 
Xu = 0.6 for our theoretical prediction which is rather low for our model and pushes it to 
the limits of its applicability. 

Finally we would like to comment on the limitations of our model and corrections to 
it. The bound state effects considered here should have received a truly non-perturbative 
treatment. We have restricted ourselves to a perturbative approximation to the problem. 
This approximation makes sense only at large enough x which we have chosen to be > 0.5. 
The contribution of more than one hard gluon exchange will be suppressed by powers of 
a s in this case. The contribution of soft gluons to the pion bound state is taken into 
account in the evolution of the distribution amplitude. The higher Fock states of the pion 
are expected to be suppressed when x is large enough [12]. The pion and the nucleon are 
not treated symmetrically in our model, namely the nucleon bound state effects are not 
taken into account since in the kinematic region we consider, x^ is always large and x u 
is always small. Gluon emission to the final state will first contribute to the evolution of 
the parton distribution functions. Hard gluon emission is suppressed by powers of a s . No 
attempt of a systematic inclusion of higher order or mass effects was made. 

The coefficient functions A, jU, and v at large x > 0.5 are very sensitive to the shape 
of the projectile's distribution amplitude 4>(z, Q 2 ), the basic hadron wavefunction which 
describes the distribution of light-cone momentum fractions in the lowest-particle num- 
ber valence Fock state. Measurements of meson form factors [12] and other exclusive 
and semiexclusive processes [16] at large momentum transfer can only provide global con- 
straints on the shape of (f>(z,Q 2 ); in contrast, the angular dependence of the lepton pair 
distributions can be used to provide local measurements of the shapes of these hadron 
wavefunctions. Detailed measurements of the angular distribution of leptons as a function 
of both x and Qt for the reactions Hp — > £ + £~X for the whole range of fixed target 
beams H = tt, K, p, p, and n will open up a new window on the structure of hadrons at 
the amplitude level. 

Our analysis shows that the broad, two-humped, distribution amplitude for the pion 
which was obtained within the context of QCD sum rules [15] can account for the main 
features of the data. In contrast, narrow momentum distributions, characteristic of weak 
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hadronic binding, predict the wrong sign for the observed azimuthal angular coefficients \x 
and v. 

Acknowledgments 

We wish to thank V. Telegdi for conversations regarding the experimental data, particularly 
the fact that the observed increase of the variable v with the dimuon transverse momentum 
is at variance with perturbative QCD calculations. We are grateful to Yu. L. Dokshitzer, 
O. Nachtmann and W. J. Stirling for useful suggestions and comments on the manuscript. 
We also benefited from discussions with H. Anlauf, V. M. Braun, T. Hyer, V. A. Khoze, 
and M. Peskin. 



7 



References 

[1] K. Freudenreich, Int. J. Mod. Phys. A 5(1990) 3643; Yu.L. Dokshitzer, D.I. 
Dyakonov and S.I. Troyan, Phys. Rep. 58 (1980) 269; Proceedings of the Work- 
shop on Drell-Yan Processes, Fermilab, Batavia, 1982 

[2] C.S. Lam and W.K. Tung, Phys. Rev. D 21 (1980) 2712 

[3] P. Chiappetta and M. Le Bellac, Z. Phys. C 32 (1986) 521 

[4] NA10 Collab. S. Falciano et al., Z. Phys. C 31 (1986) 513; NA10 Collab. M. 

Guanziroli et al., Z. Phys. C 37 (1988) 545 
[5] J.S. Conway et al., Phys. Rev. D 39 (1989) 92 
[6] J.G. Heinrich et al., Phys. Rev. D 44 (1991) 44 

[7] A. Brandenburg, O. Nachtmann and E. Mirkes, Z. Phys. C 60 (1993) 697 

[8] E.L. Berger and S.J. Brodsky, Phys. Rev. Lett. 42 (1979) 940; S.J. Brodsky, E.L. 

Berger and G.P. Lepage in Proceedings of the Workshop on Drell-Yan Processes, 

Fermilab, Batavia, 1982, p. 187 
[9] E.L. Berger, Z. Phys. C 4 (1980) 289 
[10] S.D. Drell and T.M. Yan, Phys. Rev. Lett. 25 (1970) 316 

[11] G. Altarelli, R.K. Ellis, M. Greco and G. Martinelli, Nucl. Phys.B 246 (1984) 12 
[12] G.P. Lepage and S.J. Brodsky, Phys. Rev. D 22 (1980) 2157 

[13] S.J. Brodsky, P. Hoyer, A.H. Mueller and W.K. Tang, Nucl. Phys. B 369 (1992) 519; 

K. Eshola, P. Hoyer, M. Vanttinen, R. Vogt (to be published) 
[14] S. Matsuda, Phys. Lett. B 119 (1982) 207 

[15] V.L. Chernyak and A.R. Zhitnitsky, Phys. Rep. 112 (1984) 173 
[16] A.V. Efremov and A.V. Radyushkin, Phys. Lett. B 94 (1980) 245; M.K. Chase, Nucl. 
Phys. B 167 (1980) 125; S.J. Brodsky and G.P. Lepage, Phys. Rev. D 24 (1981) 
1808; E. Braaten, Phys. Rev. D 28 (1983) 524; E. Maina and G. Farrar, Phys. Lett. 
B 206 (1988) 120; T. Hyer, Phys. Rev. D 48 (1993) 147 



8 



Figure Captions 



Fig. 1: Diagrams (a) and (b) give the leading contribution to the amplitude of reaction (4). 
Diagram (c) gives a typical (one out of four) contribution to the cross section (6). 

Fig. 2: The angular distribution coefficients A, \i and v and the Lam-Tung combination, 
1v — (1 — A), in the Gottfried-Jackson frame, versus x a for ^Q 2 T /Q 2 = 0.25. The 
dotted line corresponds to <j)(z) = 5(z — 1/2), the dashed-dotted line corresponds to 
the asymptotic <p{z) = 6z(l — z) and the dashed line shows the results for the two 
humped distribution amplitude, <p{z) = 26z(l — z)(l — 50/13 z(l — z)). The solid line 
is the result for the two-humped (f>(z) where powers of (Qt/Q 2 ) 71 ^ 2 were dropped for 



n > 3 in eqs. (13)-(15). The data points (averaged as explained in the text) are taken 
from Ref. [5]. 



Fig. 3: The same quantities as in Fig. 2 are shown, versus 




for Xu = 0.6 



and 




= 6 GeV. 



9 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9403361v2 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9403361v2 



This figure "fig3-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9403361v2 



